Delocalizing transition in one-dimensional condensates in optical lattices due to 

inhomogeneous interactions 



O 

o 

(N 



Yu. V. BludoyiQ V. A. BiazhnyiiQ and V. V. Konotopi'^.^ 

^ Centra de Fisica Teorica e Computacional, Universidade de Lisboa, 
Complexo Inter disciplinar, Avenida Professor Gama Pinto 2, Lisboa 1649-003, Portugal 
^ Departamento de Fisica, Universidade de Lisboa, 
Campo Grande, Ed. G8, Piso 6, Lisboa 1749-016, Portugal 
^ Departamento de Matemdticas, E. T. S. de Ingenieros Industriales, 
Universidad de Castilla-La Mancha 13071 Giudad Real, Spain 

It is shown that inhomogeneous nonlinear interactions in a Bose-Einstein condensate loaded in an 
optical lattice can result in delocalizing transition in one dimension, what sharply contrasts to the 
known behavior of discrete and periodic systems with homogeneous nonlinearity. The transition can 
be originated either by decreasing the amplitude of the linear periodic potential or by the change of 
the mean value of the periodic nonlinearity. The dynamics of the delocalizing transition is studied. 
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I. INTRODUCTION 

Spatial localization of the energy is a fundamental 
physical phenomenon which attracts a great deal of at- 
tention in very different physical applications and mathe- 
matical statements. When localized states are originated 
by the balance between nonlinearity and natural disper- 
sion of the medium, they are referred to as solitons. If 
localization is caused by the interplay between nonlin- 
earity and periodicity, induced either by variations of 
parameters of the medium or by discreteness, the respec- 
tive solutions are usually termed localized modes (or gap 
solitons). Localized modes are known to exist in numer- 
ous systems possessing the two main ingredients - non- 
linearity and periodicity. Discrete systems [ij, photonic 
crystals |2i] , Bose-Einstein condensates (BECs) loaded in 
optical lattices (see e.g. d, are only few examples 
where the localized modes have been discovered. 

In this context BECs in optical lattices have attracted 
a special attention, just after they have been created ex- 
perimentally 0. This is due to flexibility of optical lat- 
tices allowing one to change their parameters in time and 
in space by simple management of geometry and intensity 
of laser beams, as well as due to possibilities of experi- 
mental realization of one-, two-, and three-dimensional 
(below ID, 2D, and 3D, respectively) periodic systems. 
Naturally, a question about survivals of localized states 
subjected to variation of the lattice parameters, or in 
other words about existence of delocalizing transition 
has raised. This issue has been addressed within the 
framework of the discrete nonlinear Schrodinger equa- 
tion [gI, which models the Gross-Pitaevskii (GP) equa- 
tion in terms of the Wannier function expansion [7], and 
directly within the framework of the GP equation with a 
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periodic potential [s!] . The most important finding of the 
carried research can be formulated as the existence of the 
delocalizing transition in 2D and 3D cases and absence 
of such transition in ID systems. The phenomenon was 
explained [§] by the fact, that in ID lattices there exists 
no threshold for exciting localized states, which at small 
amplitudes become matter gap-solitons (see e.g. [1,|3,|1]). 
Such a threshold appears in 2D and 3D systems. As it 
was mentioned in [6], this explanation well corroborates 
with analogous phenomenon known for discrete systems, 
where existence of discrete envelope solitons of arbitrar- 
ily small amplitudes is possible in the ID case fl^ but in 
two dimensions requires a threshold intensity 11]. 

In this paper we show that, in a sharp contrast to the 
described situation, in the case of periodic nonlinearity 
delocalizing transition is possible in a ID system. This 
transition, resulting from the interplay between dimen- 
sionality and nonlinear periodicity, is studied in detail 
within the framework of the ID GP equation and can be 
induced by change of either the lattice potential or the 
nonlinearity. Focusing on the application of the results 
to the mean-field theory of condensed atomic gases, we 
notice that periodically varying nonlinearity appears in 
description of BEG in optical lattices which are subjected 
to optically or magnetically induced Feshbach resonance. 
Then a nonlinear lattice can be created by applying an 
optical standing wave periodically modulating the scat- 
tering length [12| , by illuminating the regions of the linear 
lattice minima by focused laser beams introducing local 
changes of the scattering length, or by a spatially peri- 
odic magnetic field which can be produced by an array of 
superconducting or ferromagnetic materials ^13] . Period- 
ical nonlinearity appears also in a natural way when one 
considers a boson-fermion mixture in an optical lattice 
provided a number of fcrmions significantly exceeds the 
number of bosons. In that case the stationary Thomas- 
Fermi distribution of the spin polarized fermions modifies 
the linear lattice and creates an effective nonlinear lattice 
for bosons due to the interspecies interactions, as this was 
discussed in detail in Refs. Il4l.[l5l. 
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The organization of the paper is as follows. In Sec. |TT] 
we introduce the model and describe the physical argu- 
ments supporting existence of the delocalizating transi- 
tion. In Sec. mil we discuss how the parameters of the 
system should be chosen and describe details of the dy- 
namics of the delocalization. Sec. II VI is devoted to a new 
mechanism of the delocalizing transition, which is based 
on the change of the nonlinearity controlling parameter 
(which is referred simply as a number of particles). The 
outcomes are summarized in Conclusion. 



II. THE MODEL 

To explain the physics of the phenomenon we start 
with the dimensionless ID GP equation 



(1) 



where U{x) is a 7r-periodic potential, U{x) = U{x + ir), 
and g{x) is the coordinate-dependent nonlinearity, which 
is also TT-periodic: g{x) ~ Q{x + 7r). In such a scaling the 
energy is measured in the units of the recoil energy Er = 
h^TT^ /{2m(f), where d is the lattice constant, and the 
spatial and temporal variables in the units d/ir and h/Er, 
respectively. Particular physical units of Q (x) depend on 
the apphcation of the model (see e.g. jl3i. 

Eq. O p ossesses stationary localized mode solu- 
tions pXTia Iit\. tp = (/>(a:)e"*'^*. Here ^J^ is & chemi- 
cal potential. It must belong either to a finite n-th gap 

{n — 1,2...): i.e. ^ G ^£ri~\fn^'^, where is either 
lower [a stands for " — " ) or upper (cr stands for " +" ) edge 
of the gap, or to the semi-infinite gap ji G ^— oo,fQ^^^ 
of the spectrum of the operator C = —d'^/dx^ + 
i.e. jO(pn ^ — En '^ip'k\ where ip'h\x) is a normalized 
Bloch state: \(pi^\x)\'^dx = 1. All branches of such 
solutions are parameterized by the chemical potential. 
They either collapse at some values of /i inside a gap or 
end up at a gap edge. The lowest branch may bifur- 
cate from the extended Bloch state. If this happens, the 
mode represents a small amplitude gap soliton governed 
by the nonhnear Schrodinger equation (see e.g. [1, Q) 

iAr = - (2M,[''^^ A^^ + Xn^\A\^A for a slow enve- 
lope amplitude A(t,£,), which is introduced through the 

representation ip w eA{T,^)ip^\x) and depends on slow 
variables ^ = ex and r — e^t. The small parameter e is 
proportional to the energy detuning towards the gap from 



the edge 6^-^^; e 



On I 



d^sk"^ /dp 



IS 



Q{x)\ipn' {x)\'^dx is the 



the effective mass and Xr 
effective nonlinearity. 

In a general situation, when the lattice amplitude de- 
creases, the width of a gap decreases as well and all gaps 
collapse at lA{x) = 0. In its turn a small gap may allow for 
the existence of only small amplitude solitons which tend 
to extended Bloch waves in the limit when the chemical 



potential approaches the gap edge. Existence of the soli- 
tons prevents delocalizing transition, because adiabatic 
decrease of the lattice amplitude leads to increase of the 
soliton width and subsequent increase of the potential 
results in recovering a localized state [1] . 

Thus, in order to obtain delocalizing transition one has 
to create a situation where small amplitude excitations 
do not exist. This is, in particular, the case of 2D and 3D 
lattices [1, H, [Til ■ In order to explore such a possibility 
in the ID model JT]) we recall that the condition for the 
modulational [1, [§] (also referred to as dynamical poj ) 

instability now reads Mn'^x^n'' < 0. For the existence of 
small amplitude solitons this condition must be satisfied 
near at least one of the gap edges. This happens, in 
particular, in models with the homogeneous nonlinearity 
{g =const) [1 i] and in discrete models with on-cite 
nonlinearity Thus, the requirements 



Mi+hi^^ > and Mt^xi-^ > 0, 



(2) 



when satisfied simultaneously, represent the conditions 
for possibility of delocalizing transition. It is to be men- 
tion here that inequalities ^ represent the constrains of 
nonexistence of stationary gap solitons, while other type 
of moving coupled-mode solitons can be excited in shal- 
low lattices when the gap width becomes small enough 
(such solitons are well known in nonlinear optics [18| 
and for different condensed atomic gases were recently 
reported in Ref. [Sill). 

The main difference of the problem at hand from the 
cases of GP equation either with a periodic potential and 
a homogeneous nonlinearity [E, Q or with periodic non- 
linearity and zero potential [l6| , is that the obtained cri- 
terion (21) strongly depends on the Bloch states and a 
specific form of the nonlinearity g{x). 



III. DELOCALIZING TRANSITION 

The signs of the effective mass at the different gap 
edges are defined: =fM^^'' < 0. Hence, in order to 
satisfy (H)) we have to check the signs of the effective 

nonlinearity X"^^ ■ We consider a stable localized state 
in a semi-infinite gap, i.e. with a solution of ([1]) with 

/i < Sq^'' ■ To be specific we explore the model with 
U{x) = -ycos(2x) and g{x) = G - cos(2x). The lat- 
tice amplitude V and the "average" nonlinearity G are 



used to control signs of the effective nonlinearities Xn 
Let us denote by G'i'^'' the values of G at which Xn''' be- 
^"^^ - 0. Since the effective nonlinearity 



comes zero, x 
is a functional of the Bloch states, it depends on lat- 
tice amplitude V, i.e. G'k^ = gIi\v). Then the curves 
G^T^ (V) on the plane (V, G) separate domains where the 
respective nonlinearities Xn^ acquire positive and nega- 
tive values (for the semi-infinite and the first bands they 
are shown in Fig. [1]). Consequently, in the semi- infinite 
gap, we have to look for the delocalizing transition in the 
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separate domains where 
Xn'iVjG) < (above and below the 



(y, G)-domain above Gq^^ In the first gap the transition 
can be expected for the parameters corresponding to the 
region above G^^'' and below g\ ' . 

Let us now suppose that G is fixed. Then, to estimate 
the lattice amplitude Vthr at which the transition occurs 
we first analyze the dependence of the number of particles 
TV = \ip{x)\'^dx on the chemical potential. For the 
lowest branch of the localized modes in the semi-infinite 
gap such dependencies are depicted in Fig. [5] for different 
V. Each of the curves has a minimum whose location 
we denote as {|I*,N^,). When the lattice potential de- 
creases, the position of the branches on the plane (/x, N) 
is changed in such a way that its minimum, = A^*(l/), 
grows. Physically this means that for localization in more 
shallow lattices higher nonlinearities are necessary. 

Let us assume that at some initial value Vini the chosen 
localized mode corresponds to the point A (see Fig. [2]), 
i.e. it has a number of atoms, which we denote Na- 
Then, one can compute the threshold value Vthr, as one 
resulting in the minimum of the numbers of particles to 
occur precisely at the same number as one of the origi- 
nal mode has, i.e. at N^,{Vthr) = Na (in the figure this is 
indicated by the point B). Designating by Vmin, the min- 
imum achieved by the lattice amplitude, one concludes 
that if Vjnin < < Vthr, the condition N^,{V) = Na can- 
not be satisfied, what in the respective dynamical prob- 
lem must result in delocalizing transition. If however 
Vmin > Vthr adiabatic increase of the lattice amplitude 
will (approximately) restore the initial form of the mode. 

The described method of finding the threshold value of 
the lattice amplitude was based on the static arguments. 
Hence the real dynamical value of Vthr can be slightly 
different, even when the change of the potential is adi- 
abatic. Therefore we provided dynamical study of the 
transition using the following variation of the potential 

amplitude: V{t) = Vrmn + iVini-Vmin)\l-2t/tf\, withtf 

being final time of simulations and V^(0) — V{tf) — Vini- 
The minimal value Vmin is achieved at tf/2. We notice, 
that considering a BEC in a nonlinear lattice created due 
to the Feshbach resonance the mentioned variation of the 
linear potential can be achieved by changing the inten- 
sity of laser beams creating the linear lattice subject to 
constant external field originating Feshbach resonance. 




FIG. 2: Number of particles vs chemical potential /i for 
the semi-infinite gap for V = 3 and V ~ 2.4. The arrow AB 
corresponds to Na = 26.25. 
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FIG. 3: (a) Dynamics of the chemical potential ij,{t) for 
Vmin = 2.4 (lower solid line) and Vmin = 2.35 (lower dashed 
line) . The upper lines show the respective changes of the up- 
per edges £q^^ of the semi-infinite gap. In (b) and (c) the 
corresponding dynamics of the density profiles are shown for 
Vmin — 2.4 and Vmin = 2.35 respectively. In (b) and (c) by 
lines we show t — tf/2 where V{t) — Vmin- Here Vini ~ 3, 
G = 0.85, tf = 1000, and N = 26.25. 



In boson-fermion mixtures the suggested situation can 
also be realized by simultaneous change of the transverse 
confinement of the gas, the latter affecting the effective 
nonlinearity together with the light beams creating the 
linear periodic potential (the details as well as examples 
with particular boson-fermion mixtures can be found in 
Ref. [ll). 

A typical result of the simulations is shown in 
Fig. [3] (a)-(c). While in panels (b) and (c) we show the 
two possible scenarios: absence of delocalization (panel 
(b)) and delocalizing transition (panel (c)), in the panel 
(a) we illustrate the evolution of the chemical potential, 
computed at each moment of time as fi = (E + Eni)/N 
where 



E = 



g{xM^ dx (3) 



is the energy of the condensate and Eni — ^ J g(x)\il;\'^dx 
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is the energy of the two-body interactions. As it is clear, 
has physical meaning only during the stage of the 
evolution until the mode is localized and follows adiabatic 
change of the potential. That is why the dashed curve 
in Fig. [3] (a) is broken at some time corresponding to the 
delocalizing transition. 

To describe the phenomenon in terms of the energy, 
we recall that for a stationary solution 

^lN = E + Eni and N = (4) 

oil 

Thus, in the semi-infinite gap, a minimum of as a 
function of /i implies maximum of the energy of a local- 
ized mode: — E{fi^) — max£'(/x). Indeed, from the 
fact that the minimum of N{fi) is achieved at /i = 
and that a given branch E is uniquely determined by N 
(and hence by /x), it readily follows that dE{^^,)/dfi — 0. 
Next, differentiating twice the first of equations in ^ 
with respect to /i twice and using the second relation, 
one obtains 

d^E_dN^^d^N 

Since in the semi-infinite gap /i < one deduces from 
([5]) that a local minimum of the number of particles with 
respect to the chemical potential corresponds to the lo- 
cal maximum of the energy of the stationary mode. We 
notice that there exists a significant difference in the en- 
ergetic properties of the modes in semi-infinite and in 
domains where the chemical potential is positive (they 
correspond to finite gaps): in the last case corre- 
sponds to the local minimum of the energy with respect 
to the chemical potential. 

The dynamical results depicted in Fig. [3] show the 
predicted delocalizing transition for Vmin < Vthr- The 
threshold value of the potential amplitude found from 
the dynamical simulations remarkably well matches the 
one found from the arguments based on the static consid- 
eration. If the minimal value of the potential Vmin > Vthr 
after increase of the potential amplitude the localized 
modes is almost identically restored. Such behavior to- 
gether with the analysis of the behavior of the chemical 
potential which approaches the gap upper edge Sq^'^ as 
V decreases [see Fig. [3] (a)], strongly support the physical 
picture of the delocalizing transition suggested above. 



magnetic field through the Feshbach resonance fl3| (what 
is possible since usually the resonant magnetic fields for 
boson-boson and boson-fermion interactions are differ- 
ent). Moreover, an interesting way of generating delo- 
calizing transition described in this section comes from 
the fact that in the quasi- ID situation, we are consid- 
ering, the quantity N, here referred to as a number of 
particles, is related to the authentic number of atoms J\f 
by the relation N = {4,asd/ {'Ka?))M (here a is the trans- 
verse linear oscillator length of the bosons). This means 
that the effective nonlinearity can be also controlled by 
the lattice period. Increase of G, i.e. the delocalizing 
transition, can be achieved by the decrease of the lattice 
constant d, which in its turn is manipulated by the angle 
between the laser beams producing the lattice [1, 0] . 

The threshold value Gthr at which the transition oc- 
curs can be found in a way similar to one we used to find 
Vthr- It is based on the requirements ^ and is schemat- 
ically depicted in Fig. [4] where we show a shift of the 
lowest branch of the solution subject to change of the 
average nonlinearity G. Since the number of particle Na 
remains constant, adiabatic increase of G corresponds to 
the shift along the line A— >B (Fig. [5]) , provided the initial 
modes is chosen in the point A. The critical value Gthr 
is found form the relation N^,{Gthr) = Na- Thus, des- 
ignating by Gmax the maximal value of G, we conclude 
that delocalizing transition occurs if Gmax > Gthr and 
does not occur otherwise. 

In order to test these arguments in the dynamical 
scheme we use the following dependence of the nonlin- 
earity on time: G(t) = Gmax + {Gini - Gmax)\'^ - 2t/t/|, 

i.e. G initially increases until t = tf /2 and then decreases 
to its initial value Gmi - The results of the simulations are 
shown in Fig. [S] In panel (a) one can observe dynamics 
of the chemical potential which at initial stages of evo- 
lution approaches the gap edge and then either returns 
to its initial value (no delocalization, solid line) or looses 
its meaning at the delocalizing transition (dashed line). 
Panels (b) and (c) show explicit behavior of the spatial 
profiles of the modes when the delocalizing transition is 
induced by the nonlinearity with Gmax > Gthr , and when 
the transition is absent, i.e. for Gmax < Gthr- As above 
the analysis of the evolution of the chemical potential 
strongly supports the suggested physical interpretation 
of the phenomenon. 



IV. NONLINEARITY INDUCED 
DELOCALIZATION 

In the case of periodically varying interactions, in- 
crease of G can also result in delocalizing transition, 
because it may change the sign of Xo^^ (see Fig. [T]). 
Such transition occurs at a constant amplitude of the 
linear lattice. As an example, we mention that such 
adiabatic change of the average nonlinearity can be 
achieved by changing of the boson-boson s-wave scatter- 
ing length fls in the boson-fermion mixture by external 



V. CONCLUSION 

In the present paper we have shown that delocalizing 
transition can be observed in one-dimensional periodic 
systems with specific spatial dependence of the nonlin- 
earity, which guarantees stability of the Bloch states at 
the upper edge of a semi-infinite gap or at the both edges 
of a finite gap. The physics of the transition is based 
on prohibition of existence of small-amplitude gap soli- 
tons bifurcating form the band edge. The delocalizing 
transition can be originated and controlled by variation 
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FIG. 4: N vs fi for the semi-infinite gap for Go ~ 0.85 and 
Go ^ 0.86. The arrow AB corresponds to iV = 26.25. 




FIG. 5: (a) Dynamics of fi for Gmax = 0.86 (lower sohd line) 
and G^aa: = 0.862 (dashed line). Now = -0.937 (upper 
solid line). Dynamics of the density profiles for Gmax ~ 0.86 
(panel b) and Gmax = 0.862 (panel c). In (b) and (c) by lines 
we show the moment t = tf/2. Here dni — 0.85, 1^ = 3, 
tf = 1000, and iV = 26.25. 



of either the lattices amplitude, or lattice period, or by 
change of the average nonlinearity, the latter achievable 
by means of magnetic of optical Feshbach resonance. In 
all the cases the phenomenon manifests itself in impos- 
sibility of restoring a localized shape of the mode if the 
change of the respective parameter goes beyond a critical 
value. The critical values of either the lattice amplitude 
or average (effective) nonlinearity can be computed with 
help of the suggested algorithm based on the analysis of 
the lowest branches of the localized modes. 



The described phenomenon is not limited to the BEC 
applications, and should be also observable on other 
physical systems. In particular, we mention that in 
the tight-binding approximation the nonlinear Gross- 
Pitaevskii equation with linear and nonlinear lattices is 
reduced to a discrete nonlinear Schrodinger equation with 
nonlocal nonlinear interactions, and thus the latter is 
probably the best candidate for observing one more real- 
ization of delocalizing transitions in one-dimension. 
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